Magneto-optical response of CdSe nanostructures 
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We present theoretical calculations of the Lande g-factors of semiconductor nanostructures using 
a time-dependent empirical tight-binding method. The eigenenergies and eigenfunctions of the 
band edge states are calculated as a function of an external magnetic field with the electromagnetic 
field incorporated into the tight-binding Hamiltonian in a gauge-invariant form. The spin-orbit 
interaction and magnetic field are treated non-perturbatively. The g-factors are extracted from the 
energy splitting of the eigenstates induced by the applied magnetic field. Both electron and hole g- 
factors are investigated for CdSe nanostructures. The size and aspect ratio dependence of g-factors 
is studied. We observe that the electron g-factors are anisotropic and find that the calculated values 
agree quantitatively with experimental data. We conclude that the two distinct g-factor values 
extracted from time resolved Faraday rotation experiments should be assigned to the anisotropic 
in plane and out of plane {qz) electron g-factors, rather than to isotropic electron and exciton 
g-factors. We find that the anisotropy in the electron g-factor depends on the aspect ratio of the 
nanocrystal. The g-factor anisotropy derived from the wurtzite structure and from the non-unity 
aspect ratio may cancel each other in some regime. We observe that hole g-factors oscillate as a 
function of size, due to size dependent mixing between the heavy hole-light hole components of the 
valence band edge states. Extension to the calculation of exciton g-factors is discussed. 
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I. INTRODUCTION 

Spin dynamics in semiconductor nanostructures have 
been studied intensively in recent years, motivated by the 
emerging field of semiconductor spintronics and quantum 
information processing.^ The most important time scale 
when implementing the quantum computer is the deco- 
herence time of the quantum degree of freedom which 
is intended to be used as the qubit. Typically the spin 
decoherence time in the bulk semiconductor material is 
extremely short. However, it is expected that the spin 
decoherence time should increase substantially in nanos- 
tructures due to the three-dimensional quantum confine- 
ment. This expectation is supported by the optical orien- 
tation experiments of Gupta et al^ where a nano-second 
spin lifetime was measured for neutral CdSe nanostruc- 
tures. This indicates that there will be plenty of time to 
perform quantum operations on the spin degree of free- 
dom in semiconductor nano structure before the coher- 
ence is lost. Consequently spins in nanostructures are 
excellent candidates for qubits. On the other hand, both 
spin based quantum computation and spintronics require 
precise control of the spin. Since the control of the spin 
dynamics in nanostructures is strongly dependent on the 
g-factors of electrons, holes, and excitons in the nanos- 
tructure, it is imperative to understand the behavior and 
magnitude of g-factors. 

Experimentally the g-factors of CdSe nanocrystals 
with wurtzite lattice structure have been measured via 
Time Resolved Faraday Rotation (TRFR)2'3 and Mag- 
netic Circular Dichroism (MCD).^ The TRFR experi- 
ment measurement reveals multiple g-factors. Two or 
four distinct g-factor values are extracted, depending on 
the size of the nanostructure which ranges from 22 A to 
80 A in diameter. The size distribution of the sample 
is about 5-15%.^ The sample from 22 A to 57 A has a 



size-dependent mean aspect ratio ranging from 1.0 to 1.3, 
with a ±0.2 variation2iiii4. It was suggested that when 
there are four distinct g-factors, these g-factors should be 
assigned to anisotropic electron and exciton g-factorsi^ 
It has been speculated that there exists a quasispheri- 
cal regime in which the g-factors become isotropic, hence 
reducing the number of distinct g-factors from four to 
two. On the other hand, MCD measurements reveal 
only a single exciton g-factor for nanocrystals with 19 
A and 25 A in diameter.^ It should be noted that the 
hole spin is also initially aligned by the optical pump- 
ing in a TRFR experiment. It has been argued that the 
fast decoherence of the hole spin makes it impossible to 
detect the hole g-factor in TRFRi^ Although it is well 
established experimentally® and theoretically^ that the 
hole spin decoherence time in the bulk semiconductor is 
extremely small, the three-dimensional quantum confine- 
ments might also alter significantly the hole spin decoher- 
ence time in the nanostructure. To the best of our knowl- 
edge there is neither experimental measurement nor the- 
oretical estimation about the hole decoherence time in 
those nanostructures. Recent time-resolved photolumi- 
ncscence on InAs/GaAs quantum dots* suggests that nei- 
ther the electron nor the hole spin relax on the lifetime 
scale of the exciton in the system but no estimation of 
the hole relaxation time could be made. This suggests 
that a hole decoherence time becomes much longer in 
nanostructures. It is thus not yet clear whether the hole 
g-factor signature should appear in a TRFR experiment. 

Theoretically, the size dependence of the electron g- 
factor in CdSe nanostructures has been calculated within 
the eight band Kane modeliMLiS, and in the tight- 
binding model. ■''3(14 In general, the effective mass approx- 
imation (EMA) type calculation is inadequate for nanos- 
tructures at small sizes ( < 30A ) because the atomic 
nature and surface effects become more prominent as the 
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size of the nanostructure decreases. Because of its atomic 
nature, the tight-binding model is ideal to study the elec- 
tronic and optical properties of nanostructures in this size 
range. 

The time-independent tight-binding approach to calcu- 
lation of g-factors in Ref 14 is based on Stone's formula^^ 
which is derived form the double second order perturba- 
tion in terms of the spin-orbit interaction and the exter- 
nal magnetic field. The results obtained from this per- 
turbative analysis show strong shape dependence. It was 
observed that a transition from anisotropic to isotropic 
g-factor tensor occurs at aspect ratio « 0.3, resulting in 
a quasispherical regime, as originally suggested by Ro- 
dina and coworkers. Since the spin-orbit interaction is 
strong in CdSe, {Xq^ ^ 0.151eV and Age = 0.320eV), it 
is desirable to include this non-perturbatively in order to 
get quantitative values. This not only provides a more 
accurate estimation of the electron g-factors for nanos- 
tructures with strong spin-orbit interaction, but will also 
enable us to treat the hole and exciton g-factors system- 
atically, in addition to the electron g-factors. It is also 
intriguing to investigate the possibility of a quasispher- 
ical regime in more a quantitative calculations for the 
electron g-factor. 

In this paper we present such non-perturbative the- 
oretical calculations of the electron and hole g-factors 
for CdSe nanostructures employing the time-dependent 
empirical tight-binding theory. Both the spin-orbit in- 
teraction and the external magnetic field are taken into 
account non-perturbatively. The g-factors are extracted 
from the magnetic field induced energy shifts of the elec- 
tron and hole eigenvalues. The size and aspect ratio de- 
pendence of the resulting g-factors is investigated. We 
observe that the electron g-factors decrease monotoni- 
cally as a function of size and are strongly anisotropic. 
The calculated values agree well with the experimental 
data from TRFR. It will be shown that the electron g- 
factors can explain the TRFR experimental data without 
invoking a exciton g-factor. The result also shows par- 
tial cancellation between the anisotropy deriving from 
the wurtzite structure and from the aspect ratio. The 
hole g-factors show very different behavior. They show 
marked oscillations as a function of the size. This is due 
to the size sensitive mixing between the heavy hole and 
light hole components. The exciton g-factor is not in- 
cluded in the current calculation. However the same cal- 
culation scheme can be easily extended to calculate the 
exciton g-factor with Coulomb interactions included non- 
perturbativelyiiSiSS 

The paper is organized as follows: In section^lwe sum- 
marize the empirical tight-binding Hamiltonian for CdSe 
nanostructures and describe how to solve the model us- 
ing a time-dependent approach. In section lTTTl we present 
our numerical results, including total density of states, 
band gap, and electron and hole g-factors. In section Hvl 
we draw several conclusions and discuss possible future 
directions. 



II. THEORY AND MODEL 

A. Tight-binding model of CdSe nanostructure 

We start from the empirical tight-binding model for 
the bulk CdSe semiconductor with an sp^s* basis. The 
parameters we use for CdSe are derived from the em- 
pirical parameters obtained by Lippens and Lannooii 
for the bulk CdSe in the zinc-blende structure, assum- 
ing nearest-neighbor interactions only. We construct the 
CdSe nanocrystals with wurtzite structure correspond- 
ing to the typical CdSe nanostructures seen in TEM 
images. The constructed structures have approximate 
but not exact C'sv symmetry. The same structures have 
been used in previous time-independent tight-binding 
studiesii^iiSiSi We remove the dangling bonds on the 
surface by shifting the energies of the corresponding hy- 
brid orbitals well above the conduction band edge. The 
spin-orbit interaction is included in the Hamiltonian. 
Spin-orbit coupling constants are assigned to both types 
of atoms, with Aq^ = 0.151eV and Agg = 0.320eV 
respectivelyi^iiS^ In order to reproduce the A-B splitting 
within the sp^s* basis, a crystal field of —40 meV is added 
to the Pz local orbitalsiSi 

B. Time-dependent approach 

Instead of diagonalizing the tight-binding Hamiltonian 
directly we employ the time-dependent approach which 
has been used previously to calculate electronic and ex- 
citonic properties of CdSe nanocrystals with zinc-blende 
structureii^ In the following we will briefiy describe the 
time-dependent tight-binding technique, indicate the dif- 
ferences compared with previous calculations, and show 
its advantage for calculation of g-factors. 

The time-dependent method depends primarily on the 
spectral decomposition for an arbitrary initial state. Let 
\En) be the complete set of eigenfunctions of the Hamil- 
tonian. Any initial state \tp{0)) can be expressed as the 
linear combination of the eigenfunctions 

\m)^Y.^n\En). (1) 

n 

The wavefunction at a later time t is 

m)) = e-'^*|V(0)) = ^5„e-^^"*|i?„). (2) 

n 

Projecting the wavefunction at time t onto the initial 
wavefunction and perform the Fourier transform one 
finds 




Thus the resulting Fourier spectrum can give us the spec- 
tral weight of the initial state in the eigenfunction basis 



3 



and the eigenenergies of the eigenstates, provided that 
the eigenstates have non-zero overlap with the initial 
state. To get the total density of states one can sum over 
the spectral decompositions obtained using each wave- 
function in a complete set as the initial state in term. 
The natural and convenient complete set to choose in 
the tight-binding framework is the direct product set of 
all local site orbitals, atomic-orbitals, and spin states. 
Then 



the expansion 



n ila ■'~ 



dte'^'{i;aamMt))- (4) 



where = \site, orbital, spin). To achieve S- 

function resolution one would need to have the infinite 
length record of the correlation function {il}{0)\ip{t)). In 
practice only a finite length T of record is available, which 
gives rise to artificial sidebands around a broadened S- 
function approximation. The finite record length is taken 
into account by multiplying the right hand size of Eq. Q 
with the normalized Hamming window function w{t)^'' 
where 



^-iHdt _ ^-i(Hi-\ hHn)dt 



e 



-iHidt 



^ ^-iHndt^-iHndt ^ 



(8) 



To implement this decomposition we first break the tight- 
binding Hamiltonian into the on-site self-energy terms, 
the local spin-orbit terms, the local Zeeman terms, and 
the hopping terms. The on-site spin-orbit interaction is 
diagonalized and exponentiated analytically in the basis 
of the tight-binding orbitals, i.e. the 6x6 matrix of 
the p-orbitals with spin. For the hopping terms we fur- 
ther use the checkerboard decompositioni2i2& to divide 
these to different independent directions. Note that in 
the zinc-blende structure there are only four fundamen- 
tal directions while in the wurtzite structure there are 
seven fundamental directions. As a result of this decom- 
position each term contributing to the short time prop- 
agator can consequently be evaluated analytically-'^^ and 
the time evolution of the state can be calculated very 
efficiently. 

The eigenfunction \En) with eigenenergy En can be 
calculated from 



w{t) = 1 -cos(2^*' 
= 0, 



if 0<t<T, 
if t > T. 



(5) 



\En) « / dte'^"*|VW), 



(9) 



The window function will reduce the sidelobes of the 
broadened 5-functions and generate a normalized peak 
height. The resulting spectrum is of the form 

(6) 

where yV„ represents the absolute spectral weight in 
eigenstate and the lineshape function ll{E — En) 

is defined by 



^i(E~E„)T+2sTT _ 



^i{E-E„)T _i ^ 

C{E-En) = ^^E~En)T "2^^^ i{E~En)T + 2sn 

(7) 

If the total wavefunction propagation time is T, the en- 
ergy resolution is ~ /S.E — tt/T. If the energy dif- 
ference between the desired eigenenergy En, and ad- 
jacent eigenenergies is larger than tt/T, the spectrum 
near energy En can be approximately represented by 
yVnC{E — En) with very high accuracy. Assuming this 
form, the value of the eigenvalue can be determined with 
accuracy much higher than tt/T. To get the most ac- 
curate value possible it is desirable to perform the time 
integration of Eq. 10 by direct integration instead of us- 
ing a discrete Fourier transform. 

In order to use the spectral method one must be able to 
calculate the time propagator e~'^* efficiently. In order 
to accomplish this we first break the time propagator 
into a series of short time propagators e~*^* = (^^-^Hdi-^N 
with t — Ndt. For the short time propagator we make 
use of the Baker-Hausdorff formula^ expansion to obtain 



provided there is non-zero overlap between the ini- 
tial wavefunction and the desired eigenfunction, i.e. 
{En\ip{0)) 7^ 0. Typically when there is no magnetic field 
(B = 0), the initial state is taken as a uniform super- 
position of local orbitals with specific angular momen- 
tum index. The resultant eigenfunctions are then used 
as the starting point to calculate the eigenfunctions and 
eigenenergies when the magnetic field is turned on. 

If there are degenerate eigenstates, the right hand side 
of Eq. (PJ will in general be some unknown linear com- 
bination of these eigenstates. However, if a set of exact 
or approximate quantum numbers which can be used to 
label the degenerate eigenstates are known in advance 
eigenfunctions corresponding to definite quantum num- 
bers can be derived by judiciously choosing an initial 
state having the same quantum numbers. Typically the 
angular momentum index is used in this work for this 
purpose. This property will be used to generate Kramers' 
doublets in our calculations. 

There are three important energy scales in this prob- 
lem. The first energy scale is the energy difference be- 
tween lowest conduction electron and higher energy con- 
duction electrons and the difference between highest va- 
lence hole and lower energy hole states. This energy scale 
is typically at the order of 100 meV or larger. The sec- 
ond energy scale is the energy difference between nearly 
degenerate hole states that correspond approximately to 
the heavy hole and light hole states in the bulk limit. 
This energy scale is size dependent and is sensitive to 
the shape of the nanocrystals. In our calculation we find 
this energy scale to be 1-100 meV. The last important 
energy scale is the magnetic field induced splitting for a 
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Kramers' doublet from which the g-factors are extracted. 
Typically this energy scale ranges from couple of hundred 
/LteV to several /xeV. 

The maximal total propagation time is about 1280000 
1/eV, resulting in an energy resolution of 2.5 /ieV. In 
our calculation this energy resolution is enough to single 
out the spectrum of band edge electron and hole states 
from other higher energy states. It is also sufficient to re- 
solve the two nearly degenerate hole states. Once a high 
resolution eigenfunction is generated, by using window 
function Eq. H5I6|I to suppress the contribution from ad- 
jacent eigenstates , the eigenenergies of band edge states 
can be determined with accuracy up to 1 /ieV. 

C. Calculation of the g-factors 

It is important to clarify the definition and the sign 
convention for the g-factors, especially when these g- 
f actors are anisotropic. When there is no external mag- 
netic field Kramers' theorem guarantees that each eigen- 
state is at least two fold degenerate. In bulk CdSe the 
heavy hole and light hole are degenerate at the F point. 
In a CdSe nanostructure it is expected that the quan- 
tum confinement will lift this degeneracy. In this work 
the g-factors will be defined with respect to the Kramers' 
doublet. For a Kramers' doublet the effective magnetic 
Hamiltonian has the form 

H^siB) = ^ibB- G -S, (10) 

where S is the effective spin operator, which is defined 

with respect to the two Kramer's states |V'*)j and G is 
the 3x3 g-factor tensor. In the Kramers' state basis, the 
effective spin operator S has the form 

where m is an integer chosen so that the real and effec- 
tive spin are approximately equal. The G tensor can be 
diagonalized in the principal axes. Let e^, i — 1,2,3 be 
the principal axes. Then an external magnetic field 

B = Biei + ^262 + ^363 (12) 

will give rise to a Zeeman splitting 

AEiB) = {E{B)-E{Q))^fiBj2g.B,. (13) 

i 

We will denote the values gi as principal g-factors. Indi- 
vidual gi can be identified via varying external magnetic 
field along each of the principal directions and calculating 
the field dependent Zeeman splitting. 

In our calculation for CdSe nanocrystal we find that 
the conduction electron is primarily s-like. In this case 
one can identify the effective spin operator with the real 
spin operator and give the g-factors a definite sign. By 



choosing the initial wavefunction to have a well defined 
real spin index, the resultant eigenfunction will approxi- 
mately have the corresponding effective spin index. 

On the other hand, we find the hole to be primarily 
p-like. Furthermore, there is a strong mixing between 
heavy hole and light hole component which is sensitive 
to the size and the shape of the nanostructures. It is 
thus difhculty to link the effective spin to the real spin, 
resulting in an ambiguity of the sign of the hole g-factors. 
In this situation it is more appropriate to express the 
Zeeman splitting by the quadratic form 

AE^iB) EE (EiB) - E{0)f ^tilY. (W) 

i 

in which the sign of principal g-factor is not well defined. 
Typically the sign convention of the g-factor in the atomic 
and bulk limit can then be used as a convention to assign 
an definite sign to the g-factors in the nanostructure. In 
both cases there usually exists a simple relation between 
the effective spin operator and real spin operator, which 
enables us to determine the corresponding sign of the g- 
factor in the nanostructure. However there is a dearth 
for experimental results of electron, hole, and exciton g- 
f actor in bulk CdSe. Hence it is difficult in this case to 
use bulk experimental results as a guide to determine the 
sign of the g-factors. 

In order to assign the sign to the hole g-factors for 
CdSe we therefore adopt the following scheme. The hole 
wavefunction will be calculated by propagating an ini- 
tial state which has definite angular momentum index 
j = '^ijz = 3/2,1/2. The sign is then determined by 
whether the magnetic field induced energy shift is posi- 
tive or negative. In the bulk limit this scheme will repro- 
duce the heavy hole and light hole value correctly. In our 
calculation we find that the electron g-factor decreases as 
the size increases. This represents qualitatively the same 
trend as seen in effective mass type calculations^ for var- 
ious semiconductor nanostructures. On the other hand 
the calculated hole g-factor for CdSe shows oscillations, 
making correlation with the atomic and bulk limits more 
difficult. 



D. Gauge invariance 

Since in this work the g-factors will be determined 
via the energy splitting of the electron and hole states 
under the external magnetic field, it is critical to cast 
the tight-binding model into a gauge invariant form. We 
use the Peierls-coupling tight-binding scheme to ensure 
the gauge invariance in our tight-binding modeliS2422i2£ 
In this scheme an electromagnetic field specified by 
the scalar potential ^{r,t) and the vector potential 
A(r, t) will modify the on-site (a, Ri|i?|a, R,) and off- 
site (a', RJ|iJ|a, Ri) tight-binding parameters via 

(a,R,|i?|a,R,) (a, R,|i/|a, R,) - $(R„ t), (15) 
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TABLE I: Size, diameters, and aspect ratio of the nanostructures. 



Number of atoms 


66 


108 


144 


237 


336 


384 


450 


561 


758 


768 


777 


1501 


\/LxLy (A) 


13.38 


13.38 


16.92 


21.85 


21.85 


25.39 


26.76 


22.85 


34.55 


27.42 


20.44 


43.01 


(A) 


11.38 


18.38 


18.38 


14.88 


21.88 


24.38 


21.88 


35.88 


35.88 


39.38 


49.88 


42.88 


Aspect Ratio 


0.85 


1.37 


1.09 


0.68 


1.00 


0.99 


0.82 


1.64 


1.04 


1.44 


2.44 


0.99 



and 

(a',R',|ff|«,R,) ^ {a' ,R'^\H\a,R,)e-'^ 

(16) 

where a straight hne should be taken for the vector poten- 
tial integral. In this calculation, the on-site intra-atomic 
dipole matrix elements {ipiialSrilipui a-') are set to zero. As 
pointed out by Foreman^, this choice ensures gauge in- 
variance but cannot correctly describe intra-atomic tran- 
sitions in the atomic limit. However, predictive calcu- 
lation of atomic transitions is not contained within em- 
pirical tight-binding treatment. Previously, optical ab- 
sorption spectra have been dealt with in tight-binding 
analysis by incorporating dipole matrix elements as extra 
fitting parameters. For CdSe nanocrystals, calculation of 
the bulk absorption spectrum has been found to be in- 
sensitive to the magnitude of the on-site dipole matrix 
elements ji^iSSiSi This provides some empirical justifica- 
tion for setting these matrix elements to zero in order to 
achieve gauge invariance. 

For gauge-invariant correction we thus only need to 
modify the the hopping constant between nearest neigh- 
bors. Since there are only seven independent hopping 
directions in a wurtzite structure the gauge phase can be 
calculated and stored before performing the time prop- 
agation. A brief summary of the gauge phase in the 
wurtzite structure is given in the appendix □ 

To estimate the contribution to the Zeeman splitting 
from the gauge phase, we have calculated the Zeeman 
splitting without gauge phase for some of the nanocrys- 
tals. We find that the gauge phase contributes 10-40% 
of the Zeeman energy. Without the gauge phase the Zee- 



man splitting increases and becomes more isotropic. 

III. RESULTS 

We investigate CdSe nanostructures having 66-1501 
atoms. This roughly corresponds to the size range of 15- 
43 A in effective diameter. In our calculation we define 
the aspect ratio to be the ratio between effective in-plane 
diameter {yJL^Xy) and out-of- plane diameter {L^). The 
aspect ratio of these nanostructures ranges from 0.68 to 
1.64. These can be divided into three different aspect 
ratio groups. The first group has aspect ratio well below 
one, ranging from 0.68 to 0.85. The second group has 
aspect ratio approximately equal to one, ranging from 
0.99 to 1.09. The third group has aspect ratio well above 
one, ranging from 1.37 to 1.64. A nanostructure with 
aspect ratio 2.44 is also studied, in order to probe the 
trends of g- factors in the quantum rod limit. In Table 
^ we summarize the in-plane diameter and out-of-plane 
diameter values of the nanostructures. If the aspect ratio 
of the nanostructures deviate from 1.0 by less than 10% 
then it is appropriate to use a single effective diameter 
[yj LxLyLz) to characterize the nanostructure. Note that 
the nanostructures used in the TRFR experiment have 
reported aspect ratios in the range of 1. 17-1. 34 f2 However 
a single effective diameter was nevertheless used to char- 
acterize the nanostructures. Furthermore, the TRFR 
sample has 5-15% size distribution and ±0.2 aspect ratio 
variation. Hence, one must be cautious when making a 
quantitative comparison between the calculated and the 
experimental results. 



To verify that the tight-binding model can reproduce 
the general features of conduction band, valence band, 
and identifiable band gap for nanostructures we have cal- 
culated the total density of states (TDOS) for smaller 
nanostructures (66-450 atoms). In Fig.^we plot the low 
resolution (« 50 meV) TDOS for a 450 atom CdSe nanos- 
tructure. It is evident from the figure that the conduction 
band edge (CBE), valence band edge (VBE), and band 
gap can be easily identified. It should be noted that the 
TDOS calculation is computationally expensive because 
one has to sum over a complete set of initial states. How- 
ever only the states at the band edges are relevant to the 
optical orientation experiment. A prior knowledge of the 
TDOS is not necessary for calculation of the band edge 



eigenstates. A reasonable initial guess of the band edge 
eigenenergy is sufficient for calculation of high resolu- 
tion band edge eigenenergies and eigenfunctions through 
an iterative procedure described below. For the smaller 
nanostructures where we have calculated the TDOS, we 
use the band edge energies identified from the TDOS data 
as initial values. For the larger nanostructures, we assign 
the initial value of band edge energies by extrapolating 
the band edge energies of the smaller nanostructures. 

To get the high resolution band edge eigenenergies 
and eigenstates we first estimate the eigenenergies as de- 
scribed above. A low resolution eigenstate is then gen- 
erated using some judiciously chosen initial state. The 
initial state is set up to have non-zero overlap with 
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TABLE II: Spectral weight of \^'^'') in the |s, (t = \) local orbital of Cd and Se. 



Number of atoms 


66 


108 


143 


237 


336 


384 


450 


561 


758 


768 


777 


1501 




0.75 


0.79 


0.81 


0.83 


0.85 


0.86 


0.86 


0.88 


0.87 


0.88 


0.89 


0.90 



to - 
Si 



450 Atom CdSe Nanostructure 



-5 -4 -3 -2 -1 1 2 3 4 5 

Energy (eV) 



FIG. 1: Total density of states for a 450 atom CdSe nanocrys- 
tal. 



the desired eigenfunction and to posses a well-defined 
value of some particular quantum number such as the z- 
component of the local total angular momentum \j,jz)- 
This low resolution eigenstate is then put through the 
spectral weight analysis described in section III Bl which 
results in a higher resolution eigenenergy. The higher 
resolution eigenenergy is then used together with the 
lower resolution eigenstate to generate a higher resolu- 
tion eigenstate. This process is iterated until the desired 
accuracy is acquired and, in the case of the hole, until the 
near degeneracy between heavy hole like and light hole 
like doublets is lifted. Once the CBE and VBE eigenener- 
gies are found, the band gap can be trivially calculated, 
from -^Gap = ^CBE ~ ^VBE- Fig. H we plot the 
high resolution results for the size dependent CBE en- 
ergy, VBE energy, and band gap. These results are all 
stable with respect to further iteration. Note that the 
VBE consists of two nearly degenerate Kramers' dou- 
blets. As the size of the nanostructure increases, these 
two doublets will converge respectively to the heavy and 
light hole doublets in bulk CdSe. 



A. Electron g-factor 

Although the CdSe nanostructures we studied here 
have only approximate but not exact C^y symmetry, we 
expect nevertheless that the principal axes are still lo- 
cated approximately along the x, y, and z directions. 



This is supported by the result of a perturbative time- 
independent tight-binding calculation of g-factors for 
these same nanocrystals^'* in which it was found that 
9x ~ Qy ^ 9z- To accurately identify the Zeeman 
splitting it is necessary to generate the Kramers' dou- 
blet which will evolve into the Zeeman eigenstates when 
we turn on the external magnetic field. (This is essen- 
tially equivalent to solving the zeroth order degenerate 
perturbation problem.) The Kramers' doublet 
for a magnetic field pointing in the +z-direction can 





















rt' . 







200 300 

Number of Atoms 



FIG. 2: (a) Band gap (b) CBE energy (c) VBE energy as a 
function of the number of atoms. Note that the VBE consists 
of two nearly degenerate levels in a CdSe nanocrystal, each 
corresponding to a perturbed Kramers' doublet, i.e. 4 states 
in total. 



be generated via setting all local orbitals of the initial 
states to have spin equal to ±i. The Kramers' doublet 
for X and y directions are then calculated as = 
± 1^-^)) and \^^y) = ± jIV'-")) re- 

spectively. The external magnetic field is limited to be 
less than 10 Tesla, which corresponds to the range of 
magnetic field in the typical experimentsSii. To make 
the connection to the CBE in the bulk material, which 
is s-like, we calculate the spectral weight of the 
state in the |s,cr — \) local orbital of Cd and Se. In 
Table ^ we summarize the size dependence of these s- 
orbital spectral weights. We find that the CBE electron 
in the nanostructure is still primarily s-like, with spectral 
weights greater than 0.75 for all sizes. The s-orbital con- 
tribution increases monotonically as the size increases. 



J 



In Fig. O we plot the magnetic field dependent spec- 
tra for the IV'^^) state with magnetic field in the x- 



direction and for the state with magnetic field in 
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=1.17 
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□A.R 
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TA.R 
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AA.R 


=1.37 


1.64 
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Energy (eV) 

FIG. 3: Zeeman shift of the l^*"*"^) component of the Kramers' 
doublet in the CBE when (a) the external field B is in the +x- 
direction and (b) the +z-direction. 



the z-direction, for a 336 atom CdSe nanocrystal. As- 
suming the spectral peaks have hneshapes of the form 
£{E — En{B)), the magnetic field dependent eigenenergy 
En{B) can then be determined with very high accuracy. 
The g- factor is then extracted by fitting En{B) as a func- 
tion of B. In Fig. ^ we plot the resulting electron g- factors 
as a function of the length parameter Lz- The data are 
grouped according to the aspect ratio of the nanostruc- 
ture. Group 1 (down triangles) has aspect ratio 0.68- 
0.85, group 2 (open squares)has aspect ratio 0.99-1.09, 
and group 3 (up triangles) has aspect ratio 1.37-1.64. 
One calculation for nanocrystal with aspect ratio 2.44 
is also included (closed circle). The extracted g- factors 
from TRFR experimental are also plotted for comparison 
(asterisks). Note that the size distribution of the sample 
in the TRFR experiments is about 5-15%, which is rep- 
resented in the figure by the horizontal error bar. The 
aspect ratio of the sample in TRFR experiment in this 
size range is about 1.17-1.34, with a ±0.2 variation. 

From the calculations we find that Qx — fjy > dz for 
all the nanostructures. As a result only two sets of data 
are shown in the figure and (jy is used to represent both 
and gy. The results show strong anisotropy between 
gz and g\\ . Both g-factors decrease monotonically as a 
function of the size of the nanostructure. The value of 
gz decreases rapidly, while the value of decreases more 
gradually. Fig.^shows that our results are in good agree- 
ment with the experimental values. It is evident from 
the figure that when there are two distinct g-factors ob- 
served for a given nanostructure, they should be identi- 
fied with the in-plane (^n) and out-of-plane (gz) g-factors 
of the electron. This assignment is very different from the 
original experimental suggestion that one of the g-factors 
should be identified with the isotropic electron g-factor 
while the other g-factor might be identified with an ex- 
citon g-factor 1^ 



10 15 20 25 30 35 40 45 50 

(Angstrom) 

FIG. 4; CdSe nanocrystal electron g-factors as a function of 
the nanocrystal length parameter L^. 



It is intriguing to look into the aspect ratio depen- 
dence of the g-factors in more detail. We observe from 
Fig. 21 that both and gz increase as the aspect ra- 
tio increases, provided that the aspect ratio is less than 
1.64. We find that gz is more sensitive to the aspect 
ratio and increases much more with this than . The 
g-factors begin to saturate between aspect ratio 1.64 and 
2.44. It is expected that if one continues to increase the 
aspect ratio then gz should begin to decrease, since it 
eventually should approach the bulk value. On the other 
hand gy should stay roughly constant after it saturates, 
provided that the in-plane cross-section is kept the same 
when one increases the aspect ratio. It should be em- 
phasized that the aspect ratio is only a simple indicator 
for the shape of the nanostructure. Two nanostructures 
with similar number of atoms and aspect ratio values 
might still have very different shape or surface struc- 
ture. From the observations above it is clear that we 
can identify a range of aspect ratios in which the effect 
of anisotropy of the wurtzite structure and that of the 
shape of the nanocrystal partially cancel each other so 
that the electron g-factors become more isotropic. In the 
cases studied here it appears that the cancellation is not 
complete. It also appears unlikely from these exact calcu- 
lation of g-factors that the cancellation will become more 
complete for large-size nanocrystals since the difference 
between gz and gy increases for larger nanocrystals hav- 
ing aspect ratio approximately unity. As a result a true 
quasi-spherical regime as predicted by an analysis per- 
turbative in spiniS*i^ in which the electron g-factors be- 
come isotropic may never reached. Under certain growth 
conditions it is possible to synthesis CdSe nanostruc- 
tures with zincblende structurCfSS. It is expected that for 
zincblende CdSe nanocrystals, if the shape of the nanos- 
tructure has high symmetry then there will be only one 
isotropic g-factori24 
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TABLE III: Spectral weight of hi state in the local orbitals \j — |,Jz)- (See text.) 



Number of Atoms 


66 


108 


144 


237 


336 


384 


450 


561 


758 


768 


777 


1501 




0.47 


0.03 


0.27 


0.82 


0.06 


0.30 


0.81 


0.06 


0.87 


0.37 


0.51 


0.88 




0.12 


0.78 


0.24 


0.02 


0.78 


0.21 


0.03 


0.79 


0.02 


0.11 


0.34 


0.03 




0.29 


0.00 


0.37 


0.04 


0.04 


0.39 


0.06 


0.05 


0.04 


0.18 


0.02 


0.04 




0.01 


0.03 


0.02 


0.00 


0.04 


0.02 


0.02 


0.02 


0.00 


0.25 


0.00 


0.00 



TABLE IV: Spectral weight of h2 state in the local orbitals \j = |, jz). (See text.) 



Number of Atoms 


66 


108 


144 


237 


336 


384 


450 


561 


758 


768 


777 


1501 


u = +i 


0.02 


0.46 


0.04 


0.46 


0.64 


0.40 


0.36 


0.53 


0.37 


0.37 


0.35 


0.04 


J. = +h 


0.06 


0.20 


0.03 


0.03 


0.10 


0.13 


0.13 


0.33 


0.17 


0.11 


0.10 


0.83 




0.08 


0.13 


0.37 


0.35 


0.09 


0.17 


0.41 


0.02 


0.37 


0.19 


0.43 


0.02 




0.68 


0.07 


0.43 


0.00 


0.08 


0.22 


0.01 


0.01 


0.01 


0.23 


0.06 


0.06 



B. Hole g- factor 

We have calculated the values of for the two nearly 
degenerate valence band edge doublets, which will be de- 
noted by hi and /i2- We define /ii(/i2) to be the highest 
(second highest) energy valence band state. To connect 
the hi and states to the heavy and light hole states 
in the bulk material we calculate their spectral weight 
in local orbitals possessing definite angular momentum 
quantum number \j = In Table-HTll (Table- Hvj) 



we summarize the spectral weight of the hi{h2) states. 
We find that for the nanocrystals in the size range we 
are interested in, the mixing between the ±| and ±i 
components is very strong. The mixing also seems to 
be sensitive to the size of the nanocrystal, without any 
clear trend emerging. As a result it becomes improper to 
rigorously identify the /ii(/i2) state with the heavy (fight) 
hole states respectively, and we therefore set m = 1 in 
Eq. l(TT)) for all hole states. 



In FigElwe plot the size dependent hole g-factors for 
those two hole doublets. The data has been regrouped 
into heavy- hole like and light- hole light states. For each 
size of the nanostructure we look at the spectral weight 
of hi and /i2 states and determine which is more heavy- 
holc-like. Similarly to the electron g-factors, the data for 
hole g-factors are also grouped by the nanocrystal aspect 
ratio. We observe that the two hole states have very 
different g-factors. Both holes show strong oscillations 
as a function of the size. It is important to clarify how 
these g-factors should approach the relevant bulk values 
when the size of the nanocrystal increases. In the bulk 
CdSe semiconductor, the valence band near the T point 
can be described by the Luttinger Hamiltonian.'^^ The 
states ||, ±|) and ||, ±i) are associated with the heavy- 
hole and light-hole respectively. If the heavy-hole and 
light-hole are really degenerate, then Eq. (|10|l is actually 
not appropriate, since the nature of the J — ^ angular 
momentum has to be taken into account. If the heavy- 
hole and light-hole are not degenerate and the mixing 
between ||,±i) and ||,±|) components are small, then 
one can use Eq. (|10|l and set to = 3 for the heavy-hole 
and m = 1 for the light- hole in Eq. Hll() . 

We speculate that the irregular mixing of local orbitals 
evident in Tables IIIII and IIVI is the cause of the hole g- 
factor oscillation in nanocrystal size. From the table we 
observe that the hi state becomes more and more heavy- 



hole like for nanocrystal with more than 450 atoms, pro- 
vided that the aspect ratio is close to one. It is expected 
that as the size of the nanocrystal increases, the two 
hole states will eventually converge to the heavy-hole and 
light-hole respectively. One must then be careful when 
comparing the hole g-factors values calculated here with 
the bulk heavy-hole g-factor since the latter is usually 
defined with m = 3 in Eq. ((TT)l >j^ 



IV. SUMMARY AND DISCUSSION 

We have calculated the g-factors of the conduction 
band edge electrons and valence band edge holes in 
CdSe nanostructurcs using the time-dependent tight- 
binding method. This allows an exact, non-perturbative 
analysis, with extremely high resolution of the Zeeman 
shifts. We observe that the electron g-factors are strongly 
anisotropic, with gx = gy > 9z for all nanocrystal 
sizes. The size dependence of the anisotropic electron 
g-factors agree quantitatively with the values extracted 
from TRFR experimental data. This leads to the con- 
clusion that when there are two distinct g-factors in the 
TRFR experiments, they should be identified with the 
electron in-plane ((7||)and out-of-plane ((7z)g-factors re- 
spectively. This is very different from the original exper- 
imental suggestion that one of the g-factors is an isotropic 
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FIG. 5: CdSe nanocrystal hole g-factors as a function of the 
nanocrystal length parameter Lz- 

electron g-factor while the other one is an exciton g- 
factoriS We have investigated the aspect ratio depen- 
dence of the electron g-factors. We find that in general 
the g-factors initially increase as a function of the aspect 
ratio, and that increases more than g^^ . The increase 
of g-factors saturates around aspect ratio 1.62-2.44. It is 
expected that gz will then begin to decrease again until 
it reaches the bulk value while will stay roughly the 
same, if the aspect ratio continue to increase. The aspect 
ratio dependence allows us to identify a regime where 
the anisotropy derived from the wurtzite structure and 
that derived from the shape of the nanocrystal cancel 
each other partially, resulting a more isotropic regime. 
However a full cancellation is never reached in our cal- 
culation, unlike the previous observation from a pertur- 
bative analysis^ It also appears unlikely that in larger 
nanocrystals the cancellation will become complete, since 
the difference between gz and for an unit aspect ratio 
nanocrystal increases as a function of the size. 

We find that the valence band edge consists of two 
nearly degenerate Kramers' doublets, i.e. only a small 
perturbation from the bulk states. The hole g-factors 
for these states show oscillations as a function of the 
size. We speculate this is due to the strong size sensi- 
tive heavy/light hole mixing of the two hole states. 

One possible extension of the current calculation 
scheme is the evaluation of the exciton g-factors in these 
nanostructures. The g-factor of a uncorrelated electron- 
hole pair can be approximated by gx = ge — Qh- How- 
ever the Coulomb interaction gives rise to excitonic ef- 
fects and this simple picture then breaks down. The 
time-dependent tight-binding method has been success- 
fully applied to calculate excitonic properties in CdSe 



nanocrystaiia and it appears feasible to extend the cur- 
rent scheme to now calculate the corresponding exciton 
g-factors. Since the exciton fine structure splitting in 
CdSe nanocrystal is of the order 1-10 meV^^ the energy 
resolution obtained here (« 1 /iV)is enough to resolve the 
exciton fine structure. 
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APPENDIX 

In the wurtzite structre there are 7 independent elec- 
tron hopping directions. In this calculation these 7 hop- 
ping directions are denoted by ddir^ dir = 1 • • • 7, and are 
assigned to be the following vectors: 

di = f (0,0,3), 

d2 = f (2V2,0,-1), ds = f (-2V2,0,-1), 
da = f (-^/2,^/6,-l), rfg = f (A -V6, -1), 
d4 = f (-V2,-V6,-l), d; = f (V2,V6,-1). 

(A.l) 

Here ao = 2.625A is the lattice constant. This convention 
enables us to calculate the gauge-dependent quantities 
explicitly. 

For a fixed external magnetic field B = {Bj., By, Bz) 
we assign the vector potential to be yl — ^B x f. We 
define a magnetic-field-dependent gauge phase 

g niia+ddir 

(t>{RaJd^r) = - / A-dl, (A.2) 

where Ra represents the position vector of an anion. The 
gauge phase for the cation can be easily calculated by 
taking the appropriate complex conjugation of the phase 
of the corresponding anion. With the line of integra- 
tion to be a straight line, the integral can be calculated 
analytically. Using this notation the gauge-dependent 
short-time propagation of the electron hopping in some 
particular direction becomes 

g-V-,„,.,„dt (^Ai\^i{Ra)) + Am\^m{Rc))) (A.3) 

= (Ai cos{Vi„Jt) ~ iAra sin{Vi,ndt)e+'^''^^^'''"^) l*;(i?a)) 
+ (A^cos{Vi^dt) ~ azsin(l/^™di)e""^^'^" |*m(i?c)), 
where Vim is the hopping constant in zero magnetic field. 
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